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ABSTRACT 

We consider quantum correlation functions of the antiferromagnetic spin-| Heisen- 
berg XXZ spin chain in a magnetic field. We show that for a magnetic field close 
to the critical field h c (for the critical magnetic field the ground state is ferromag- 
netic) certain correlation functions can be expressed in terms of the solution of the 
Painleve V transcendent. This establishes a relation between solutions of Painleve 
differential equations and quantum correlation functions in models of interacting 
fermions. Painleve transcendents were known to describe correlation functions in 
models with free fermionic spectra. 



PACS-numbers: 75.10.Jm 02.30. Jr 



1 Introduction 



In this letter we continue our investigation of zero-temperature correlation functions of the 
XXZ Heisenberg model in the critical regime — 1 < A < 1 in an external magnetic field. The 
XXZ hamiltonian is given by 

« = E + + A - !) - K > ( L1 ) 

where the sum is over all integers j, L is the length of the lattice, <r Q are Pauli matrices and h is 
an external magnetic field. For later convenience we define A = cos(2r/), where | < 77 < 7r. The 
free fermionic point in this notation is 77 = The model ( [Lip can be solved by means of the 
Bethe Ansatz, which yields a description of the spectrum and eigenstates (with iV down spins 
and L — N up spins) in terms of the roots of the following set of coupled algebraic equations 

!§ 

/ sinh(A 3 - -iy) \ L _ sinh(A fc - Xj + lirj) - = l N - ^ 

Ysinh(Aj + M7) / sinh(Afc — Aj — 2zr?) ' 

For the case A > — 1 it was proved by C.N. Yang and CP. Yang in |J that the ground state is 
characterized by a set of real Xj subject to the equations (|1.2|) . In the thermodynamic limit the 
ground state is described by means of an integral equation for the density of spectral parameters 
P(A) 



27rp(A) — / d\x K(X, fi) p{jj) 

J-A 



sin(2r^) 



sinh(A — irf) sinh(A + irj) 



K(/i, A) = — M4r] [ f (1.3) 

smh(/i — A + 2irj) smh(/i — A — zirj) 

Here the integration boundary A is a function of the magnetic field h. The physical picture 
of the ground state is that of a filled Fermi sea with boundaries ±A. The dressed energy of a 
particle in the sea is given by the solution of the integral equation 

e(A) - - / dii K(X, n) e(n) =2h- — - \ -. S ■ (1.4) 

Z7T J —A smh(A — ir]) smh(A + ir\) 

The condition that the dressed energy vanishes at the Fermi edge e(±A) = determines A as a 
function of the magnetic field h. It was shown in that the ground state of (|1.1| ) for |A| < 1 
is partially magnetized for magnetic fields h < h c = 4cos 2 (?]). For larger magnetic fields h > h c 
the ground state is the saturated ferromagnetic state. For h — >• the integration boundary A 
tends to 00, whereas for h —>■ h c A — >• 0. Below we only consider the region h < h c and in 
particular the limiting case h — > h c . 

The subject of this letter is the generating functional of correlation functions 

G(m) = (exp(aQ 1 (m))> = (0| exp(af^ ar a +)\0) = (0| exp(a^ -^-)|0) , (1.5) 



i=i 3=1 2 



where |0) is the ground state. 
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Various correlation functions can be obtained from (exp(aQi(m))) , (e.g., the Ferromagnetic 
String Formation Probability |10 which corresponds to setting a = — oo) or 



exp(aQi(m))|0) + 1 - 4 / dX p(X) , (1.6) 

a=0 J-A 

where A is the lattice laplacian acting on a function f(j) defined on the lattice as A/(j) = 
f(j) + f(j — 2) — 2f(j — 1) and p(X) is defined in ( |1.3|) . In what follows we will consider 
G(m) in the limit h — > /i c . As is shown below G(m) is connected to the solution of a Painleve 
V transcendent. The connection of Painleve transcendents and integrable models with free- 
fermionic spectra is well established [l], [3], [2j. However, we want to emphasize that in the 
present case we are dealing with a theory of interacting fermions, so that the connection is 
novel. 



(0KoJ|0) = 2A(0| 



cr 
dc7 2 



2 Determinant Representation 



In a recent paper || we used the approach invented in || [7| (for a detailed exposition of this 
method see ||) to derive the following representation of G(m) in terms of the determinant of 
a Fredholm integral operator. 



(6| det (l + v)\0) 
(0|exp(ag 1 (m))|0) = V ' ,} > . 

det (i _ -K) 

Here K and V are Fredholm integral operators with kernels given by (|1.3j) and 

V( ^ A n = sin(2?7) f 1 + e^(A)e 2 (^) 

2tt sinh(A — fJ,)\ sinh(A — + 2irf) sinh(A — fj, — 2in) 

+ exp(a + ^) - ^(A))( -^ 1 + . f}''"'^' 

\smh(A — /i — 2«?7) smh(A — /i + 2z- 

where 



(2.1) 



e 2 (A) 



sinh(A + ^)\ m eV2(A) ^ ei(A) = f sinh(A-^ ^ l(A) 



(2-2) 



(2.3) 



,sinh(A — zr/) 7 1 ' Vsinh(A + irj) . 

The quantities <Pj(X) are bosonic quantum fields || (they are called "dual quantum fields") 
defined by 



¥>a(A) = p tt (A) + g (A) , (0| 9a (A) = = p o (A)|0) ,(6|0) = 1 , o=1...4 



[ft(^)>Pa(A)] 



/l 1 0\ 

10 1 

111 

Vi o i ij 



ab 



^10 1^ 
110 
10 11 

Vo i i ij 



ln(h(tM,\)), 



(2.4) 



ab 



where h(\, u) = &inh ^ n ^ 2r ' ? ^ . Note that the dual fields have the important property that they 
commute 



[<p a (X),<Pb{fi)] = • 



(2.5) 
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3 The limit of strong magnetic field 



For large magnetic fields h — > h c ,h < h c the integration boundary A tends to zero according to 

1, 



A = -\tm(ri)Wh c -h + 0(h c -h). (3.1) 

This fact can be used to essentially simplify the above representation as was first done for a 
different correlation function in [K|. We first expand the kernel V(A,/i) for small |A|, < A. 
We define y — m\ cot(r/)| as a shorthand notation. We then observe the following simplification 
of the commutation relations between the "momenta" p a (X) and "coordinates" g a (A) 



/ 








i 


-1\ 










-l 


1 




-i 


i 








V 


i 


-i 





o 1 



i cot(2r/)(/j - A) + C(A 2 ) . (3.2) 



ab 

This allows us to reduce the number of dual fields from four to two via the identification 
<^ 2 (A) = —ifi(X) and ^(A) = —<p 3 (\). Furthermore the r.h.s. of the commutators ( |3.2| ) is 
a linear function of A — /i so that we can choose a representation such that v 9 a(A) are linear 
functions in A 

<A»(A) = ¥a + <p' a X , Pa(A)=p a +plA, q a (A) = q a + q' a \ , a = 1, 3 (3.3) 

where [p a> g 6 ] = 0, [pj,, q' b ] = 0, and 

[ft,Ps] = -i cot (277) = -[^3,Pi] , [g 3 ,i?i] = -i cot (277) = -[gi,j/ 3 ] . (3.4) 

As v ? i(A) appears in ( |2.2| ) only in the combination ipi(\) — <fi((J>) the quantity <p\ drops out. 
This in turn implies that p' 3 and q' 3 commute with all remaining operators and thus will not 
contribute to the expectation value w.r.t. (0| and |0). Therefore we can drop them everywhere. 
In the next step we perform a similarity transformation with exp(i\(y + |^)), which leaves 
the determinant of the Fredholm integral operator invariant but brings the kernel to a more 
symmetric form, in which the dual fields now only enter via the expressions a := a — 2if 3 and 
x := y+ %(pi. In what follows it is crucial that the quantum operators in the dual bosonic Fock 
space x and a are commuting objects and no problems with operator orderings occur before 
we evaluate the expectation value w.r.t. (0| and |0). 

Putting now everything together we arrive at the following simplified representation valid 
in the limit h — > h c 

(exp(aQ 1 (m))) = (0| det (l + V ) |0) (l + ^- cot(2r/) + £>(A 2 )) , (3.5) 



where the kernel of Vq is given by 



V (X, n) = (exp(S) - 1) Sin( , ( ^ ^ . (3.6) 

Here we use the following notations for the dual fields 

a = a + a q + a p , x = y + x q + x p , (0\x q = (0\a q = , x p \0) = a p \0) = , 
[x 9 ,S p ] = [S 9 ,x p ] = — cot(2r/) , [a,x] — . (3.7) 

All other commutators vanish. We would like to emphasize that a and x commute, which is 
important for the further analysis. 
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4 Connection with Painleve V 

Let us define a new variable £ = Ax and consider the object 

a (t) =t^ln(det(l + Vj,))- 
at 



(4.1) 



In || it is shown that o"o obeys a Painleve V differential equation in the case where a and 
x are real numbers. In our case a and x are quantum operators, but due to the fact that they 
are commuting we still can follow through the derivation of ||. We thus find that o"o obeys the 
following nonlinear differential equation 



£ 



cPaojty 
dt 2 



dt 



dt 



(4.2) 



which is identified (see |§) as a r - function form of the Painleve V equation. Rewriting 
in terms of the function yo(t) defined through 



u 2 (t) 

a {t) = -Aitu{t) + — y-(y (t) - l) 2 , u(t) 



4tty (t) - £gg 

Hy (t) - 1) 2 



(4.3) 



one obtains (see again ||) standard form of the Painleve V differential equation for the function 

w(*) = y (|) 



d 2 u{t) _ fdu(t)\ 2 3u(t)-l 2u(t)(u(t)-t 
dt 2 ~ I dt J 2u(t)(u(t) - 1) + u(t) - 1 



1) 2m;(£) ldivjt) 

t t dt 1 ' ' 



The large £ asymptotics of the solution of the above equations are known (see [fj, |TT|, [TB|, [□[]) and 
can be used to extract the large distance asymptotics of our Fredholm determinant. Combining 
the results of [0 (asymptotics for cr ) and [[II], |I^] (constant term) we obtain that 



ln(det(l + y )) 



2 a 2 , . . . „. , ,a,, a 3 

-Axa + — In 4Ax + 2 In <? — + — — 
it 2tt z 2n Sn^xA 

a 2 . 5a 4 

■ cos(46») 



1 



x- 



32tt 2 A 2 



128tt 4 A 2 



+ 0( 



A 3 y 3 



(4.5) 



where the distance y must be large (the product of y and the small parameter A should go to 
infinity) and where 



9\y) 



,(l+7)^ 2 



e « 



n=l 



AO 



exp ( v~ — — I dt 

— J 



2 ia 
4Ax + -Sln(4Ax) - 4argr( — ' 

7T 2"/T' 



+ iVt) + - iVt) 



(4.6) 



where ^(z) = 4-\aT(z) is the digamma function and 7 is Euler's constant. Formula (|4.6|) for 
the phase 9 of the corresponding solution of the Painleve V equation (|4.4j) was also obtained in 



12| . In order to obtain the large-distance asymptotics of the correlation function we still have 



to evaluate the expectation value in the dual bosonic Fock space. It is easiest to evaluate the 



quantities {Q\{m) k ) = ( 



gh 



3 aQi(m) 



o=0 



) directly. Expanding 



ln(4Az) 



C3 



—2 

e27r- 



• ln(4Ay) 



1 + 



Ot, OC ' -p I X q 

2tt 2 y 



+ 0(y 



1 + 



a 



+ 0(A^y~ 2 ) , 



8ir 3 Ay 

we obtain the leading terms in the asymptotic decompositions as is shown in the Appendix 



(4.7) 



(Qi(m)) 
((Qi(m)) 2 } 



lead 



lead 



m\ cot rj\ 1H cot(2r/) + 0(A 2 



71 



71 



< 2A , ,s 2 

— m| cot 7/1) + 



ln(4Am| cot 77 1) 



7T 



7P 



(l + 8Acot(2,) +0(A2)) 



+(l+ 7 )(— + 



1 8Acot(2?7) 6Acot(2^) , 



7T 

+ 0(A 2 ) 



(4.8) 



7T 7r° 71"° 

Analogous formulas can be derived for (Qi(m)) k for k > 2. Using ( |1.6|) we are now in a position 
to determine the large-distance asymptotics of the (cr^crl) correlation functions. By acting with 
(twice) the lattice laplacian on (|4.8|) we obtain 



m 2 Tc 2 



1 



3 A cot (2??) 



7T 



C(A 2 



(4.9) 



This does not yet include the contribution from the cos-term in (|4.5|) , which is very small 
as far as ((Qi(m)) 2 ) is concerned but becomes important upon differentiation. The leading 
contribution can be obtained by the methods explained in the appendix and is found to be 



7T" 



+ 0(A) cos (4Am| cot rj\{l + 0(A))) 



m" 



where 



8Acot(2r?) ^ /l9 . 
= 2 + ^ + 0(A 2 ). 

7T 



(4.10) 



(4.11) 



Our final result for the first three terms of the asymptotics of the correlation function is thus 



IA\ cot (77) I 



I A cot (2f])' 



TV 



m 2 7T 2 



+ 



7T 



7T 2 m e 



cos(4Am| cot 77 1) + 



(4.12) 



where the errors in the A-expansion are given above. This agrees with the result obtained by 
means of finite-size corrections and Conformal Field Theory. 
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5 Summary and Conclusion 



In this letter we have established a connection between the generating functional of correlation 
functions G{m) ( |1.5|) (and thus the correlator (cr^cr^)) of the spin | Heisenberg XXZ model in a 
magnetic field close to critical and the Painleve V differential equation. Painleve transcendents 
were known to describe correlation functions for models with free-fermionic spectra |], |2|. For 
generic magnetic fields in the XXZ spin chain the general approach of || should be followed: 
the determinant representation ( |2.1|) , ( |2.2| ) should be used to embed the quantum correlation 
function in an integrable system of integro-difference equations and one then should solve the 
associated Riemann-Hilbert problem. 
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A Appendix 

In this appendix we discuss how to evaluate the expectation value with respect to the dual 
quantum fields and explicitly evaluate the quantities ((Qi(m)) k ) . According to 



((QM)) k ) 



det (l + Vo] 10) [ 1 + — cot(2r/) + 0(A 

a=0 V ' V VT 



(A.l) 



Using ( |4.5| ) for the asymptotics of the logarithm of the determinant and expanding according 
to (}4.71 ) we see that in order to get the leading asymptotics we need to evaluate the derivatives 
of expectation values of the form 



(0|e 



2 A ^ e ^l n( 4A y) 



1 + 



01. p I CC q 
2^2 y 



1 + 



a 



8n 3 Ay 



|0). 



(A.2) 



Using the commutation relations ( |3.7|) we see that we can bring all terms into the form (we use 
that (0\x q = 0) 

(0\x™e^ A ™F(a)\0) , (A.3) 

where F(a) is only a function of a and contains no x p 's or ir q 's. The central identities we will 
use in order to evaluate the expectation values are 



z m = (0\x™e« Axa F(a)\0) 



(A.4) 
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where k = 1 — ^-cot(2rj). The identities are established via induction. The induction start 
m = is proved as follows. Expanding the exponential and using that (0\x q = 0, S p |0) = 
and [a p , a q ] = (to move all S p 's to the right) we obtain 



00 1 /2A\" 

*>=(0|E-7( — ) (y + x p ) n (a + a q rF(a)\0) . (A.5) 



n=0 U - V 71 



Expanding 



77 ' 

(6|(i/ + %)> + a,)" = £ ' ^, (oi(i/ + %Y*> n - k (A.6) 

fc=0 «!(n - K)\ 

and then using the commutation relations 

(6|[/(2 P ),S g fc ] = (0|(cot(2r ? )) fc / (fc) (2 P ) , (A.7) 
where is the k'th derivative of the function /, we arrive at 



n=0 fc=0 



(A.i 



We now use the integral representation 



1 Ire* 

dt—rr, (A.9) 



(n-ife)! 2m J t n - k+1 

where the integration contour is a small circle around the origin (and we integrate in the 
mathematically positive direction) in order to be able to perform the ^-summation (which is 
of the form of a binomial sum) 



1 r p* 00 /2A\" ~ n 
-« = ^/*jE(— ) (0|[oot(2»7) + T (y + ^)]»f(a)|0) 

(A.10) 



2m „ „ n=0 



The n-summation can be performed using (1 — z) 1 = J2h^=o z ■ Finally we perform the t- 
integration formally using the identity 

— fdt 6 —- = e°(«> , (A.ll) 

where 0(a) is an operator depending on a. Here we need to keep in mind that we are interested 
in evaluating derivatives with respect to a at a = 0. This yields the result ( |A.4j ) for m — 0. 
The induction step goes as follows. We rewrite z m as 

z m = (Ol^-^e^^c^lO) , (A.12) 

where we used that x p |0) = 0. Evaluating the commutator using 

x p /(a)|0) = cot(2r / )/'(a)|0) , (A.13) 
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where /' is the derivative of / and collecting terms we obtain 



z. 



„, - -(() ./• 

K 



m— 1 
P 



cot(2r])e^ xa F'(a) + — cot(2r])ye^ xa F(a) 

7T 



10) 



(A.14) 



Using the induction assumption and again ( |A.13| ) then yields the desired result ( |A.4|) . Let 
us now demonstrate how to evaluate the leading contributions to (|A.2|) for (Qi(m)) and 
{(Qi{m)) 2 ) . They are given by 



((Qi(m)) 1 ) 



lead 



d l 
da 1 



(0\e^ Axa e& ln(4Ay) (g(—)) 2 \0)[l + —cot(2r l ) + O(A: 2 ) 



2A 

TV 



(A.15) 



-~2 



We apply (|A.4j) with m = and F(a) = ^^(g^)) 2 , then perform the differentiations 
with respect to a and set a to zero, and finally use (|A.13 ) to evaluate the expectation value 
using the fact that (0\a q = 0, a p \0) = 0. The function g has the properties that g(0) = 1, 
g'{0) = and ^"(O) =2(1 + 7), where 7 is Euler's constant, which leads to the result 



(Qi(m)) 



2A 



lead 



lead 



7T 



6A 



— m\ cot r] I H cot (27/) + 0(A 



TV 



, 2A , ,x 2 

( — m|cotr/|J + 

7T 



ln(4Am| cot ry|) + 1 + 7 



7T- 2 



8Acot(2r7) ^ /a9nn 
x(l + ^JL + 0{k 2 )). 

IT 



(A.16) 



The contributions of the subleading terms can be taken into account in an analogous way, 
which leads to the results ( |4.8| ). We note that contributions from further subleading terms are 
of higher order in A and y. 



8 



References 

[1] E. Barouch, B.M. McCoy, T.T. Wu, Phys. Rev. Lett. 31 (1973) 1409, 
C.A. Tracy, B.M. McCoy, Phys. Rev Lett 31, 1500 (1973), 
T.T. Wu , B.M. McCoy, C.A. Tracy, E. Barouch, Phys. RevB13, 316 (1976), 
B.M. McCoy, C.A. Tracy, T.T. Wu, J. Math. Phys. 18 (1977) 1050, 
M. Jimbo, T. Miwa, Proc. Japn. Acad. 56A (1980) 405, 
B.M. McCoy and J.H.H. Perk, Nucl. Phys. B285 (1987) 279. 

[2] B.M. McCoy, S. Tang, Physica 20D, 187 (1986). 

[3] M. Jimbo, T. Miwa, Y. Mori, M. Sato, Physica ID, 80 (1980) 

[4] R. Orbach, Phys. Rev 112, 309 (1958). 

[5] C.N. Yang, CP Yang, Phys. Rev 150, 321 (1966), Phys. Rev 150, 327 (1966). 
[6] V. E. Korepin, Comm. Math. Phys. 113, 177 (1987). 

[7] A. R. Its, A. G. Izergin, V. E. Korepin, and N. A. Slavnov, Int. J. Mod. Phys. B 4, 1003 
(1990). 

[8] V. E. Korepin, A. G. Izergin, and N. M. Bogoliubov, Quantum Inverse Scattering Method, 
Correlation Functions and Algebraic Bethe Ansatz, Cambridge University Press, 1993. 

[9] F. H. L. EBler, H. Frahm, A. G. Izergin, and V. E. Korepin, Commun. Math. Phys. 174, 
191 (1995). 

[10] F. H. L. EBler, H. Frahm, A. R. Its, and V. E. Korepin, Nucl. Phys. B 446 [FS], 448 
(1995). 

[11] E. Baser, H. Widom, Jour. Fund. Analysis 50, 387 (1983). 

[12] B. Suleimanov, Lecture Notes in Mathematics vol. 230, 1191 Springer Verlag 1986. 
[13] A.M. Budylin, V.S. Buslaev, Algebra and Analysis 7, 74 (1995). 



9 



